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J£ ' Abstract 

Q : 

We consider a family of pseudo differential operators {A + <z Q A Q / 2 ; a <G [0, 1]} on R d that 
evolves continuously from A to A + A"/ 2 , where d > 1 and a £ (0, 2). It gives rise to a family of 
Levy processes {X a , a G [0, 1]}, where X a is the sum of a Brownian motion and an independent 
symmetric a-stable process with weight a. Using a recently obtained uniform boundary Harnack 
principle with explicit decay rate, we establish sharp bounds for the Green function of the process 
X a killed upon exiting a bounded C 1,1 open set D C M. d . As a consequence, we identify the 
Martin boundary of D with respect to X a with its Euclidean boundary. Finally, sharp Green 
function estimates are derived for certain Levy processes which can be obtained as perturbations 
of X a . 
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1 Introduction 

Discontinuous Markov processes have been intensively studied in recent years due to their impor- 
tance both in theory and applications. In contrast to the diffusion case, the infinitesimal generator 
of a discontinuous Markov process in M d is a non-local (or integro-differential) operator. Most of the 
recent studies have concentrated on discontinuous Markov processes (and corresponding integro- 
differential operators) that do not have a diffusion component. See [5j[9] and the references therein 
for a summary of some of these recent results from the probability literature. For recent progress 
in PDE literature, we refer the readers to [HI El IE]- 

However, Markov processes with both diffusion and jump components are needed in many 
situations, like in finance and control theory. See, for example, [2H [23 ES] ■ On the other hand, 
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the fact that such a process X has both diffusion and jump components is also the source of many 
technical difficulties in investigating the potential theory of X. The main difficulty in studying 
X stems from the fact that it runs on two different scales: on the small scale the diffusion part 
dominates, while on the large scale the jumps take over. Another difficulty is encountered at the 
exit of X from an open set: for diffusions, the exit is through the boundary, while for the pure 
jump processes, typically the exit happens by jumping out from the open set. For the process X, 
both cases will occur which makes the process X much more difficult to study. 

Despite these difficulties, significant progress has been made in the last few years in under- 
standing the potential theory of discontinuous Markov processes with both diffusion and jump 
components. Green function estimates (for the whole space) and the Harnack inequality for some 
processes with both diffusion and jump components were established in [30\ I32| . The parabolic 
Harnack inequality and heat kernel estimates were studied in [M] for the independent sum of a 
Brownian motion and a symmetric stable process, and in [13] for much more general diffusions 
with jumps. Moreover, an a priori Holder estimate is established in [13] for bounded parabolic 
functions. Very recently, the boundary Harnack principle for some one-dimensional Levy processes 
with both diffusion and jump components was studied in [26] , where sharp estimates on Green 
functions of bounded open sets of K were also established. Most recently, a boundary Harnack 
principle with explicit decay rate for nonnegative harmonic functions of the independent sum of a 
Brownian motion and a symmetric stable process in C 1,1 open sets in M. d was obtained in |10| . 

The main goal of this paper is to use the boundary Harnack principle obtained in [10] to establish 
sharp Green function estimates in C 1 ' 1 open sets for the Levy processes that are independent sums 
of Brownian motions and symmetric stable processes. 

Let us now fix the notation and state the main result of this paper. Throughout this paper, we 
assume that d > 1 is an integer and a G (0, 2). Let X° = (X®, t > 0) be a Brownian motion in M. d 
with generator A = X^f=i anc ^ ^ Y = (Yti t>Qi) he an independent (rotationally) symmetric 
a-stable process in W d . For a > we define the process X a = {Xf, t > 0) by X® = X® + aY t , 
called the independent sum of a Brownian motion and a symmetric a-stable process with weight a. 
This process, although very specific, serves as a test case for general Markov processes with both 
diffusion and jump components. 

Let D be a C 1 ' 1 open set in R rf , let X a,D be the process X a killed upon exiting D and let 
Gp(x,y) denote the Green function of X a,D (for precise definitions see Section 2). Our main goal 
is to establish sharp two-sided estimates for G a D (x,y). Let Sd{%) denote the Euclidean distance 
between the point x G D and the boundary dD. The main result of this paper is the following 
theorem. Here and in the sequel, for a, b G R, a A b := min{a, b} and a V b := max{a, b} Define for 
d > 3 and a > 0, 




when x, y are in the same component of D 



when x, y are in different components of D; 



for d = 2 and a > 




when x, y are in the same component of D, 
when x, y are in different components of D; 



(1.1) 
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and for d = 1 and a > 0, 




(5 D (x)5 D (y)) L ^A^ 
a a ((5 D (x)5 D (y)) 1/2 A 



Theorem 1.1 Zei M > 0. Suppose that D is a bounded C 1,1 open set in M. d . There exists 
C = C(D, M, a) > 1 such that for all x,y G D and all a G (0, M] 



Note that the above estimates are uniform in a G (0, M]. In case d = 1, a (non-uniform) estimate 
is covered by [26j . Letting a J, in (jl.3p recovers the Green function estimates for Brownian motion 
killed upon exiting D; for the latter, see [TU p. 182] for d = 2 and [37] for d > 3, respectively. 

The rest of the paper is organized as follows. Section 2 gives preliminary and background 
materials. Theorem 11.11 is proved in Sections 3, 4 and 5. The proof of the theorem in the case 
d > 3 is by now quite standard. Once the interior estimates are established, the full estimates in 
connected C 1 ' 1 open sets follow from the boundary Harnack principle by the method developed by 
Bogdan [1] and Hansen |20j . However this method is not applicable when d < 2 since Brownian 
motion is recurrent in this case. When d = 2, we use a capacitary argument and some recent results 
on subordinate killed Brownian motions, which are given in Section 4. The case d = 1 is dealt with 
in Section 5, where we follow the arguments of [26J. In Section 6, using the boundary Harnack 
principle and our Green function estimates, we show that both the Martin and the minimal Martin 
boundary of the process X a,D can be identified with the Euclidean boundary when D is a bounded 
C > open set. In the last section, we extend our results on X a to symmetric Levy processes that 
can be obtained from X a through certain perturbations. In particular, for every m > 0, we obtain 
sharp Green function estimates of A + m — (m 2 ^ a — A) a in any bounded C 1 ' 1 open set with zero 
exterior condition. The process corresponding to A + m — (m 2//a — A) a is a Levy process that is 
the independent sum of a Brownian motion and a relativistic a-stable process with mass m. 

Throughout this paper, we use the capital letters C\, C2, • ■ ■ to denote constants in the statement 
of results, and their labeling will be fixed. The lowercase constants c±, c?, ■ ■ ■ will denote generic 
constants used in proofs, whose exact values are not important and can change from one appearance 
to another. The labeling of the constants c\, C2, • • • starts anew in every proof. The dependence of 
the constant c on the dimension d and a G (0, 2) may not be mentioned explicitly. The constant 
a G (0, 2) will be fixed throughout this paper. We will use ":=" to denote a definition, which is 
read as "is defined to be". B(x,r) denotes the open ball in M. d centered at x with radius r > 0. 
Recall that for any x G D, 5d(x) denotes the distance between x and dD, and for a, b G M, 
a A b := min{a, b} and a V b := max{a, b}. We will use d to denote a cemetery point and for every 
function /, we extend its definition to d by setting f(d) = 0. Lebesgue measure in M d will be 
denoted by dx. For a Borel set A C M. d , we also use \A\ to denote its Lebesgue measure. 



C- 1 g a D (x,y)<G a D (x,y)<Cg a D (x,y). 



(1.3) 
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2 Preliminaries 



A (rotationally) symmetric a-stable process Y = (Y t ,t > 0,P x ,x G R d ) in R d is a Levy process 
with the characteristic exponent i.e., 



E, 



ii\Y t -Y ) 



for every x£K d and £ G 



The infinitesimal generator of Y is the fractional Laplacian A Q / 2 , which is a prototype of non-local 
operators. The fractional Laplacian can be written in the form 

A a / 2 u(x) = lim / (u{y) - u(x)) f& dy, 

£ i° J{ym d :\y-x\>e} \x-y\ a+a 

whereA(d,-a) := a2 Q - 1 vr- a! / 2 r(^)r(l- f)" 1 . Here T is the Gamma function defined by T( A) := 
J^t x - l e~ l dt for every A > 0. 

Suppose A is a Brownian motion in M. d with generator A = Yli=i J^?) an d Y is a symmetric 
a-stable process in M. d . Assume that A and Y are independent. For any a > 0, we define the 
process X a = (A t a , t > 0) by Af := A 4 ° + aY t . As already mentioned, the process X a is called the 
independent sum of the Brownian motion X° and the symmetric a-stable process Y with weight 
a. It is a Levy process with the characteristic exponent $ a (£) = |£| 2 + a a |£| a , £ G R , and its 
infinitesimal generator is A + a a A a / 2 . The process X a has a jointly continuous transition density 
that will be denoted by p a (t, x, y). From the Chung-Fuchs criterion (see [H Theorem 1.17]), it easily 
follows that, when a > 0, X a is transient if and only if a < d, while it is well known that A is 
transient if and only if d > 3. 

There is another representation of the process X a which will be useful in Sections 3, 4 and 5. 
It can be obtained by subordinating A with an independent subordinator T" := t + a 2 Tt where 
T = (T t , t > 0) is an a/2-stable subordinator, i.e., the processes (A", t > 0) and (Aj. a , t > 0) 
have the same distribution. Note that the Laplace exponent of T a is </> a (A) = A + a a \ a / 2 . Let 
M a / 2 (t) := E^oC -1 )"^ 2 / 1 ^ 1 + W 2 )- ft follows by a straightforward integration that 

OO 1 

e- xt M 1 _ a/2 (a 2a « 2 - a k)dt = —— 
o 0°(A) 

which shows that the potential density u a of the subordinator T a is given by 

u a (t)=M 1 ^ a/2 (a 2a ^- a h). (2.1) 

Since, for any a > 0, (fi a is a complete Bernstein function, we know that (see, for instance, |30j) 
u a (-) is a completely monotone function. In particular, u a (-) is a decreasing function. Since 
u a (t) = u l (a 2a l( 2 ~°^t), we see that a \— > ii a (i) is a decreasing function. Moreover, since the drift of 
T a is equal to 1, we have that u a (0+) = 1 and so 

u a (t) < 1 for t > 0. (2.2) 

The Levy measure of X a has a density with respect to the Lebesgue measure given by 

J a (x,y) := f(y - x) := f(\y - x\) = a a A(d, -a)\x - y\^ d+a \ (2.3) 
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which is called the Levy intensity of X a . It determines a Levy system for X a , which describes 
the jumps of the process X a : For any non-negative measurable function / on R + x M. d x M d with 
f(s,x,x) = for all s > and x G M. d , and stopping time T (with respect to the filtration of X a ), 



E, 



Y,f{s,X a s-,Xs) 



E, 



T 



f{s,X a s ,y)J a {X a s ,y)dy)ds 



(2.4) 



(see, for example, Proof of Lemma 4.7] and [121 Appendix A]). 

The quadratic form (£ a ,.F) associated with the generator A + a Q A Q / 2 of X a is given by 



T = W 1,2 (R d ) := <j u G L 2 (M d ; dx) : ^- G L 2 (R d ; dx) for every 1 < i < d 



■2/iad 



du 



2/inid. 



and for u, v G J 7 , 



£ a (u,v) 



Vu(x) ■ Vv(x) dx + 



1 



dxi 



(u(x) - u{y))(v(x) - v(y)) ^ dl "J " dxdy. 

i \X — y\ a + a 



In probability theory, the quadratic form (£ a , W 1,2 (K d )) is called the Dirichlet form of X a . Let 
£f(u, u) := £ a (u,u) + f Rd u(x) 2 dx. Note that for every a > 0, there is a positive constant c = 
c(a, d, a) > 1 so that 



(\Vu{x)\ 2 + u{xf) dx < £i(u,u) < c / (|Vtt( x)\ 2 + u(x) 2 } dx 

Jm. d 



for u G W 



Thus the processes X a , a > 0, share the same family of sets having zero capacity. 

For any open set D C M. d , rfy := inf{t > : Xf ^ D} denotes the first exit time from D by 
X a . We denote by X a ' D the subprocess of X a killed upon leaving D. The infinitesimal generator 
of X a ' D is (A + a a A a ^ 2 )\o- It is known (see [13]) that X a ' D has a continuous transition density 
pf)(t,x,y) with respect to the Lebesgue measure. 

Definition 2.1 A real- valued function u defined on R d is said to be 

(1) harmonic in D C M d with respect to X a if for every open set B whose closure is a compact 
subset of D, 



E, 



\u{X% 



< oo and u(x) = E r 



u{X%) 



for every x G B; 



(2.5) 



(2) regular harmonic in D C M d with respect to X a if it is harmonic in D with respect to X a and 



u(x)=E a: u(X% 



for every x G D; 



(3) harmonic for X a ' if it is harmonic for X a in D and vanishes outside D; 

(4) superharmonic in D C M rf with respect to X a if for every open set B whose closure is a compact 

subset of D, 



E, 



|u(X?a)| < oo and u(x) >E X u(X* a ) 



for every x € B. 



(2.6) 
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It follows from [13] that every harmonic function u in D with respect to X a is continuous in D 
and/ Rd m y )|(lAH-( d+Q ))dy<oo. 

Using the parabolic Harnack inequality from [13} Theorem 6.7] and a scaling argument, the 
following uniform Harnack principle was established in [10] , 



Proposition 2.2 (Uniform Harnack principle) Suppose that M > 0. There exists a constant 
C\ = Ci(a, M) > such that for any r E (0,1], a E [0, M], xo E M. d and any function u which is 
nonnegative in R rf and harmonic in B(xo,r) with respect to X a we have 

u(x) < C\u{y) for all x,y E B(xo,r /2). 

We recall that an open set D in M. d with d > 2 is said to be C 1,1 if there exist a localization 
radius R > and a constant A > such that for every Q E dD, there exist a C ^-function <p = 
<f> Q : R^ 1 -> R satisfying 0(0) = 0, V0(O) = (0, . . . ,0), ||V0||oo < A, |V0(x) - V</>(y)\ < A\x - y\, 
and an orthonormal coordinate system CSq: y = (2/1, ... , J/d-i, J/d) = : (y, with its origin at Q 
such that 

B(Q, R) n D = {y = (y, y d ) E 5(0, i?) in CSq : y d > 4>{y)}. 

The pair (R, A) is called the characteristics of the C ' open set -D. Note that a C ' open set may 
be disconnected. Observe that the distance between any two distinct connected open components 
of D is at least R. By a C 1, open set in M we mean an open set which can be written as the union 
of disjoint intervals so that the minimum of the lengths of all these intervals is positive and the 
minimum of the distances between these intervals is positive. Note that a C ,l open set may be 
unbounded. It is well known that any C 1 ' 1 open set D satisfies the uniform exterior ball condition: 
There exists R > such that for every z E dD, there is a ball B z of radius R such that B z C (D) c 
and dB z n dD = {z}. Without loss of generality, throughout this paper, we assume that the 
characteristics (R, A) of a C ' open set satisfies R = R. 

Observe that for any C 1, open set with C > characteristics (R, A), there exists a constant 
k E (0, 1/2], which depends only on (R, A), such that for each Q E dD and r E (0, R), D n -B(Q, r) 
contains a ball B(A T (Q), nr) of radius «r. In the rest of paper, whenever we deal with C ' open 
sets, the constants A, R and k will have the meaning described above. 

Let Q E dD. We will say that a function u : M — ► R vanishes continuously on D c n B(Q, r) if 
it = on D c n -B(Q, /■) and it is continuous at every point of dD n -B(Q, r). 

The following theorem is the main result of [lOj . 



Theorem 2.3 (Uniform Boundary Harnack Principle) Suppose that M > 0. For any C 1, 

open set D inM. d with the characteristics (R, A), there exists a positive constant C2 = C2(ce, d, A, R, M) 
such that for all a E [0, M], r £ (0,1!], Q £ dD and any nonnegative function u in R rf that is har- 
monic in D n B(Q,r) with respect to X a and vanishes continuously on D c n B(Q,r), we have 

44 < c 2 ^44 ever y x > y e D n r / 2 )- ( 2 - 7 ) 

A subset I? of R rf is said to be Greenian for X a if X a,D is transient. A Greenian set for X° 
will be simply called Greenian. As mentioned in the second paragraph of Section 2, when d > 2 
and a > 0, any non-empty open set D C R rf is Greenian for A" a ; and any non-empty open set in 
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R rf is Greenian when d > 3. An open set D C M. 2 is Greenian if and only if D c is non-polar (or 
equivalently, has positive capacity). In particular, every bounded open set in R 2 is Greenian. 

For any a > and any Greenian open subset D of R d for X a , we use Gp(x,y) to denote the 
Green function of X a,D , i.e., 

poo 

G a D (x,y) := / p a D (t,x,y)dt (2.8) 



o 



where P£,(t, x, y) is the continuous transition density of X a,D with respect to the Lebesgue measure. 
The function GJ,(-, •) is finite off the diagonal. It follows immediately from (|2.8p that Gf)(x,y) is a 
positive continuous symmetric function off the diagonal of D x D such that for any Borel measurable 
function / > 0, 



D 



f(Xf)ds 



D 



G a D (x,y)f(y)dy. 



We set G a D equal to zero outside D x D. The function Gf)(x, y) is also called the Green function of 
X a in D. For any iGO, Gp( • , x) is superharmonic in D with respect to X a , harmonic in D \ {x} 
with respect to X a and regular harmonic in D \ B(x, e) with respect to X a for every s > 0. 

Recall that a point z on the boundary 3D of an open set D is said to be a regular boundary 
point for X a if P^t^ = 0) = 1. An open set D is said to be regular for X a if every point in dD is 
a regular boundary point for X a . It is easy to check that every C 1, open set D is regular for X a 
for all a > and using the argument in the last paragraph of the proof of [16\ Theorem 2.4], we 
conclude that for any bounded C 1 ' 1 open set D, G a D { ■ , z) vanishes continuously on dD for every 
z £D. 

Now, as a corollary of the uniform boundary Harnack principle and the fact that, for any 
bounded C 1 ' 1 open set D, Gp(-,z) vanishes continuously on dD for every z £ D, we have the 
following proposition. 

Proposition 2.4 Suppose that M > 0. For any bounded C 1,1 open set D in M. d with the character- 
istics (R,A), there exists a positive constant C3 = C^(a, d, A, R, M) > 1 such that for all Q £ dD, 
r G (0, R) and a E (0, M] we have 

G a D (x, Zl ) G a D (x,z 2 ) , , 

G%(y, Zl ) ~ 3 G a D (y,z 2 y { ' ' 



when x,y G D\ B(Q,r) and zi,z 2 G D n B(Q,r/2). 

The following scaling property will be used below: If (X^ ,D ,t > 0) is the subprocess in D of 
the independent sum of a Brownian motion and a symmetric a-stable process in M. d with weight a, 
then (AA^I^,t > 0) is the subprocess in XD of the independent sum of a Brownian motion and a 
symmetric a-stable process in ~R d with weight aA^ a_2 ^ a . So for any A > 0, we have 

V a \D (*> y) = A"VMA- 2 t, A" 1 ^ A" 1 ?/) for t > and x, y G AD. (2.10) 

By integrating the above equation with respect to t, we get that when D is Greenian for X a , 

G a D (x,y) = X d - 2 G a x i a - 2)/a (Xx,Xy) for x,y G D. (2.11) 

In particular, for d = 2, we have 

G2,(a;,y) = G^ Q " 2)/<1 (Ax,Ay) forx,j/GD. (2.12) 
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3 Higher Dimensional Case: d > 3 

In this section we assume that d > 3. We will use G a (x, y) = G a (y — x) = G^ d (x, y) to denote the 
Green function of X a in M. d . 

Recall that u a is the potential density of the subordinator = t + a Tf given in (|2,ip . The 
Green function G a of X a is also given by the following formula ([30]) 

roc 

G a (x)= / (Airt)- d / 2 e-^ 2/< - u) u a (t)dt, x£R d . (3.1) 
Jo 

Using this formula, we can easily see that G a is radially decreasing and continuous in M d \ {0}. 
Lemma 3.1 Suppose that M > 0. For all a £ [0, M], we have 

G M (x) < G a (x) < r^a=2, for all x £ 

Proof. We have seen that for all t > 0, u a (t) < 1, and the function o i— ► u a (t) is decreasing on 
M + , cf. (12. 2p and the text preceding it. The desired inequalities follow immediately from these 
properties and (13. □ 



Lemma 3.2 Suppose that M > 0. There exist R\ > and Li > 2 smc/i /or aZZ a € [0, M] 

G a (x) > 2G a (L lX ), for all \x\ < R 1 . (3.2) 

Proof. By |30|, Theorem 3.1], there exists c\ = c\(a,d,M) > such that 

lim G M {x)\x\ d - 2 = ci. (3.3) 

x|— >0 

Let Li = (2/ci) 2 /( d ~ 2 ) V2 and take < Si < ci(l - L~( d ~ 2 )/ 2 ). Using ([53]) . we can choose a positive 
constant i?i > such that 

( c i - 5 i)r4z9 ^ gM ( x ) when < ^i- (3-4) 



x 



Thus, by Lemma [3TTI for every \x\ < Ri 



G a (x)>G AI (x)>(ci-5i)-^>L d /-^ > 2G a (L lX ). 



□ 



The next proposition gives the interior estimates for G a L 



D- 

Proposition 3.3 Suppose that M > 0. For any bounded and connected C 1,1 open set D in M. d 
there exists a positive constant C4 such that for every a £ [0, M] 

G a D (x, y) < C 4 , ^—^ for all x,y £ D (3.5) 

\x — y\ a z 

and 

G a D (x,y) > Cl 1 _ 1 when 2\x - y\ < S D {x) A S D (y). (3.6) 

\x y\ 
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Proof. Since G a D (x,y) < G a (x,y), (|3.5p is an immediate consequence of Lemma [3TT1 So we only 
need to show (|3.6p . Without loss of generality we assume that 5d(u) < Sd(x)- 

Recall that L\ > 2 and R\ are the constants from Lemma 13.21 By Lemmas I3.1H3.21 and (|3.4p . 
we have 



G a (x,y)-G a (L lX ,L iy ) > ^G a (x,y)> l -G M (x,y) > cy \ 
for some positive constant C\. 

Case 1: L\\x — y\ < 5o{y)- We consider three subcases separately: 
(a) 5d{v) < Ri- Note that, since L\\x — y\ < <5d(?/), 

\X% -y\> 5 D (y) > L^x-yl 

B(y.S D {y)) 

Thus by the fact that G a (-) is radially decreasing and (|3.7p . 



when \x—y\ < R± (3.7) 



G a D (x,y) > G a B 



B{y,6 D (y)) 



(x,y) = G a (x,y)-E x G a {X a T , 



B(y,S D (y)) 



,v) 



> G a (x,y)-G a (L 1 x,L 1 y) > c v 



\x-y 



d-2 ' 



(b) doiy) > R\ and L\\x — y\ < R\. In this case, \X% — y\ > R\ > L\\x — y\ and, again 

B(y,Ri) 

by the fact that G a ( ■ ) is radially decreasing and (|3.7p . 



> G a (x,y)-G a (Lix,L l2/ ) > c r 



\x-y 



d-2 - 



(c) Sd(u) > Ri and Li|x — y\ > R\. In this case, we have 8d{x) > $d{v) > — y| > -Ri- 
Choose a point u> G dB(y, R\/(2L\)). Then from the argument in (b), we get 

G a D (w,y) > ci 



d-2 ' 



(Ri/{2L 1 )) 

Since D is a bounded and connected C ' open set and \x — io| < \x — y\ + \y — w\ < Sd(v)/ ' L\ + 
R\/(2L\) < 3(5d(2/)/(2Li), by Proposition 12.21 and a chain argument, we have 



G a D (x,y) > c 2 G a D (w,y) > c 3 



1 



> c 3 2 d - 2 (2L 1 /R 1 



(Ri/(2L 1 )) 



d-2 



\2(d-2) 



1 



F — y 



d-2 ' 



Case 2i 2\x - y\ < 5 D {y) < L x \x - y\. Take x G dB(y,8 D {y)/{L l + 1)). Then 



I a; - y| < ^£>(y) < Li|x - y\ = - L ] - 5 D (y) < S D (y) A 8 D (x ). 

Z Li\ -f- 1 

Since D is a bounded and connected C 1 ' 1 open set and |xo — x| < |xo — j/| + — x\ < ( £ +1 + 2)^0(2/)) 
by Proposition 12.21 a chain argument and the argument in the first case, there are constants 
Ci = Ci(D, a, L\, M) > 0, i = 4, 5, 6, such that 



G a D {x,y) > c 4 G a D (x ,y) > c 5 - 



Fo - y 



d-2 



> C6, 



\x-y 



d-2 " 
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This completes the proof of the proposition. □ 

Suppose that D is a bounded and connected C 1 ' 1 open set in M. d with characteristics (R, A) and 
corresponding k. Fix zq G D with kR < 5d{zq) < R, and let E\ := kR/2A. For x,y £ D, define 
r(x, y) := 5d(x) V <$o(y) V |x — y\ and 

y) := jz G £> : (5d(^) > ^r(a3, y), |t - z| V |y - z| < 5r(x, y) j 

if r(x,y) < £1, and B(x,y) := {zq} otherwise. 
Put C 5 := C 4 2 d - 2 (^(,z )- d + 2 . Then by (|33D . 



G^(-,z )<C 5 onD\B(z ,5 D (z )/2). 

Now we define 

<f(x) := G a D (x,z )AC 5 . 

Note that if 5d(z) < 6ei, then |z — zo| > <5d(zo) ~~ 6ei > 8d(zo)/2 since 6ei < 6d(zq)/^, and 
therefore y a (.z) = G a D {z,zo). 

Using the uniform Harnack principle (Proposition I2.2| ) and Proposition 12.41 the following form 
of Green function estimates follows from [201 Theorem 2.4]. 

Theorem 3.4 Suppose that M > 0. For any bounded and connected C 1,1 open set D in M d , there 
exists Cq = Cq(D, M, a) > such that for all x,y in D and all a G (0, M] 

6 g a (^) 2 ^^'W - 6 y a (A) 2 1 y| 

where A G £>(x, y). 

Suppose -D is a bounded and connected C > open set. For all x,y G D, we let and Q y be 
points on dD such that 5_d(x) = \x — Q x \ and 5rj(y) = \y — Q y \ respectively. It is easy to check 
that if r(x, y) < S\ 

A r (x,y)(Qx), Ar(x,y)(Qy) G B(x,y). (3.8) 

(Recall that, for any Q G dD, A r (Q) is a point such that B(A r (Q), nr) C D n B(Q,r).) Indeed, 
by the definition of ^4 r (x,j/) r(x,s/)(Qx)) > Kr(x,y) > nr(x,y)/2. Moreover, 

\ x ~ A(x,y)(Qx)\ < \x - Q x \ + \Qx ~ A(x,y)(Qx)\ < 5 D (x) +r(x,y) < 2r(x,y) 

and \y — A-(x,y) (Qx)\ <\y~x\ + \x-A r 

(x,y)(Qx)\ < 3r(x,y). This verifies the claim Q3.8p. 
Recall the fact that g a (z) = G a D (z,zo) if Sd(z) < 6ei. By Theorem 12.31 and the fact that 
CQ 1 r(x,y) < 5D(A r ( x ^(Q y )) < cor(x,y) for some cq > 1, there exists c\ > 1 such that for every 
a G (0, M] and all x,y £ D with <5d(x) < 6ei and 5d{v) < 6ei, 

1 r(x,y) ~ g a {A r ( Xi y)(Q x )) G a D {A r{x ^ y) (Q x ), z ) ~ l r(x,y) 

and 

-i Sp(y) < 9 a (y) = G a D (y,z ) ^ ei fo(y) ^ 
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Proof of Theorem 11.11 when d > 3: First we assume that D is connected. 

Combining inequalities (|3.9p and (|3.10p with Proposition 13.31 Theorem 13.41 and the fact that 



SD(x)S D (y) < / a 5 D (x)5 D (y) ^ < 9 5 D (x)5 D (y) 



(r(x,y))< 



\x -y\ 



4 (r(x,y)¥ 



(3.11) 



(see [3]), we get the inequalities (j 1 . 3 f) . 

Next we assume that D is not connected. Let (R, A) be the C ' characteristics of D. Note that 
D has only finitely many components and the distance between any two distinct components of D 
is at least R > 0. Assume first that x and y are in two distinct components of D. Let D(x) be the 
component of D that contains x. Then by the strong Markov property and the Levy system (|2.4p 
of X a , we have 



cr D (x,y)=n x 



D\D(x) 



f(\X«-z\)G a D (z,y)dz)ds 



Consequently, 



f(dmm(D))E x [r a D{x) } / < G%(x, y) < f(R) E x [r a D{x) ] / G a D (y,z)dz. 

JD\D{x) JD\D(x) 

(3.12) 

Applying the two-sided estimates (|1.3p established in the first part of this proof to D(x), we get 



<D{x) 



c\^D{x){x) = ci5 D {x) 



c\ 1 Sd(x) = c x 1 5 D{x) (x) < E x 
for some c\ = ci(D, M,a) > 1. Clearly, using (|3.13p . 

G a D {y,z)dz > f G a D(y) (y,z)dz = E y [rf ){y) ] > c 2 5 D (y). 



(3.13) 



D\D(x) 



On the other hand, it follows from (|3.5p that sup zeD ae ^ AI j E 2 [tq] < C3 < 00. Moreover by (|3,13p 
and the Levy system (|2.4p of X a , 



( G a D (V,Z)dZ < ^[Ti] =Ey 

JD\D(x) 



r D(y) 



Ex T « [r a D ] 

D(y) 



< C4 L 5 D (y)+E y 



f D{V) ( I f(\X s -z\)E z [r a D ]dz)ds 

JO \JD\D(y) J 



< c 4 6 D (y) + c 5 j M (R)E y [Tory)] < c 6 S D (y) 



We conclude from the last three displays, (|3.12p and the form of j a given in (|2.3|) that there is a 
constant C7 = Cf(D, M, a) > 1 such that for every a G (0, M], 



c 7 a a 8 D (x)5 D {y) < G D (x,y) < c 7 a a 5 D (x)5 D {y). 



(3.14) 



Since for x and y in different components of D, R < \x — y\ < diam(L>), we have established (|1.3|) . 
Now we assume that x, y are in the same component U of D. Applying (|1.3|) to U we get 



G a D (x,y)>G a u (x,y)>c 8 . 



1 

x-y 



d-2 



1 A ^t/(g)^(j/) 



\x - y\ 



C 81 



1 

\x-y 



d-2 



1 A S D (x)5 D (y) ' 



\x - y\ 
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For the upper bound, we use the strong Markov property, the Levy system (|2.4p . and (|3. 13j) — f)3. 14j) 
to get 



G a D (x,y) = GUx,y)+E x G a D {X%,y) 



< Cg- 

< Cg- 

< Cg- 



1 



x — y\ d 2 
1 

x — y\ d ~ 2 
1 

x — y\ d ~ 2 



1A^M| +El 



1 A 
1 A 



\x — y\ 2 

5p{x)5 D (y) \ t , M 

\x - y\ 2 
5p{x)8 D (y) 

\x - y\ 2 



D\U 



f(\X a s -z\)G D {z,y)dz\ds 



+ j M (R)E x [r^] [ G a D (y,z)dz 
Jd\u 

+ c w 5D{x)SD(y) / 5 D (z)dz. 
Jd\u 



(3.15) 



Since the boundedness of D implies 

5 D {x)S D {y) < en 



x-y\ 



d-2 



1 A 



we have from (|3.15p 



G a D (x,y) <c 12 - 



x - y\ c 



1 A 



5 D {x)5 D (y) 
\x - y\ 2 

Sp{x)5 D (y) 
\x - y\ 2 



□ 



Define 



a(x,y,z,w) := < 





when 


X 


y, z, w are in the same component of D, 


—a 


when 


x 


y G D(x),z £ D(x) and w G D(z), 


a 


when 


X 


weD(x),#({y,z}nD(x)) = l, 


a 


when 


X 


y,z,w are all in different components of D 


2a 


when 


X 


w G D(x), {y, z] n D(x) = 0. 



(3.16) 



The next theorem will be used in Section [71 

Theorem 3.5 (Generalized 3G theorem) Suppose that M > 0. For any bounded C ' open set 
D in M. d , there exists a constant Cg = Cg(D, a, M) such that for all x,y, z,w G D and a G (0, M\, 

G a D (x,y)G a D (z,w) 



G a D (x,w) 
< C$a(x,y,z,w) 



x — w\ A \y — z\ \ ( \x — w\ A \y — z\ 

1 VI — ■ V 1 



\z — w\ 



\x - y\ 

Proof. Recall that r(x,y) = 5d{x) V 6d{u) V \x — y\ and let 

1 S D (x)5 D (y) 



x — w 



d-2 



x — y\ 2 \z — W 



d-2 ' 



(3.17) 



9D(x,y) :-- 



\x - y\ d 2 (r(x,y)) 2 
By Theorem 11.11 for the case d > 3 and (|3.11|) , 



and H(x,y,z,w) :- 



\x — w 



d-2 



\x — y\ d 2 \z — w\ 



G a D {x,y)G a D {z,w) g D (x,y)g D (z,w) 
< cia[x,y,z,w)- 



G a D (x,w) 



9d{x,w) 

, jD{y)^D{z)r{x,w) 2 

cia(x, y, z, w) ^—H(x, y, z, w) 

r(x, y) z r{z, wy 



(3.18) 
(3.19) 
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1. If \x-w\< 8 D (x) A6 D (w), g D (x,w) > \x - w\~ d+2 . Thus, by (I3TT8D 

G a D (x,y)G° D (z,w) 



G a D (x,w) 



< aa(x, y, z, w)H(x, y, z, w). 



2. Note that if y = z, since r(x, w) < 2r(x, y) + 2r(y, w), 

8D(y)S D {y)r{x, w) 2 < 8 / 5 D (y) 2 + 5 D (y) 2 



r(x,y) 2 r(y,w) 2 \r(y,w) 2 r(x,y) 2 

Thus 

9D{x,y)g D {y,w 



9d(x,w) 



<8H(x,y,y,w). (3.20) 



Now consider the case \y — z\ < 5o(y) A 5d(z). In this case gi)(y,z) > \y — z\ d+a . Thus, 
using (|3.20|) . we obtain that 

9D{x,y)gD{z,w) g D (x,y)g D (y,z) g D (x,z)g D (z,w) 1 



< 64 



9d(x,w) go(x,z) go{x,w) go{y,z) 

\x — z\ d ~ 2 \x — w\ d ~ 2 1 



x — y\ d 2 \y — z\ d 2 \x — z\ d 2 \z — w\ d 2 gr>(y,z) 



64 \x-w\ d 2 1_ ^ 21 ^ 

\x — y\ d ~ 2 \y — z\ d ~ 2 \z — w\ d ~ 2 goiy, z) 



Thus, by (f3A8l) and (|3T211 . we have 

G a D (x,y)G a D (z,w) 
G a D (x,w) 



< c 2 a(x, y, z, w)H(x, y, z, w). 



3. Now we assume that \x — w\ > 5d(x) A 5d(w) and \y — z\ > <5o(y) A 5d(z). Since 5d(x) V 
8d(w) < 5d(x) A 5d(w) + \x — w\, using the assumption Sd{%) A 5d(w) < \x — w\, we obtain 
r(x,w) < 2\x — w\. Similarly, r(y,z) < 2\y — z\. By (|3.19|) . we only need to show that 



8 D (y)5 D (z)r(x,w) 



2 



f\x — w\A\y — z\ \ / \x — w\ A \y — z\ \ .„ 

V1 X v1 )- (3 ' 22) 



r(x, y) 2 r(z, w 

Since r(rr, to) < 2r(x, y) + 2r(y, io) < 2r(x, y) + 4r(y, z) + 4r(z, u;), we have 

5 D (y)8 D (z)r(x,w) 2 < c f 5 D {y)5 D {z) 5 D (y)5 D (z) | d D (y)8 D (z)r(y, z)' 
r(x,y) 2 r(z,w) 2 ~ \ r(z,w) 2 r(x,y) 2 r(x,y) 2 r(z,w) 2 

< c ( + + r (^' z ) 2 

yr(z, u;) r(x,y) r(x,y)r(z,w) 

( r{y,z) r(y,z) r(y,z) 2 

< C4 —. r H t r + 



v r(z,io) r(x,y) r(x,y)r(z,w) 
which is, by [221 Lemma 3.15], less than or equal to 



2c 4 



r(y,z) v A / r(y,z) v 1 



13 



On the other hand, clearly 



S D (y)^D(z)r(x,w) 2 5 D (y)5 D (z) r(x,wf 



r(x, y) 2 r(z, w) 2 r(x, y)r(z, w) r(x, y) 2 r(z, w 



2 



( r(x, w) \ ( r(x, w) 

\r(x,y) ) \r(z,w) 

Thus 

5D(y)5D(z)r(x,w) 2 ( r(y,z) Ar(x,w) \ ( r(y, z) A r(x, to) 



< c 5 ^ VI ^ V 1 

r(x, y) z r(z, wy \ r(x,y) J \ r{z,w) 

Now applying the fact that r(x,w) < 2\x — w\, r(y,z) < 2\y — z\, r(x,y) > \x — y\ and 
r(z,w) > \z — w\, we arrive at (I3.22|) . 

We have proved the theorem. □ 

Note that, since we consider disconnected open sets too, we can not apply \22\ Theorem 1.1] 
directly to get the generalized 3G theorem. 

Taking y = z in (|3.17j) . we get the classical 3G estimates, that is, 

G%(x,z)G%(z,w) ^ _ , , \x-w\ d ~ 2 ^ ^ /,,2« ^ \x-w\ d ~ 2 

D ' a D \ ' < C 8 a(x,z,z,w)- l - r ^ < C 8 M 2Q V1 i l ~ r ^ 

G]j{X,W) \x — z\ a z \z — w\ a z \x — z\ a z \z — w\ a z 



4 Two Dimensional Case 

In this section we assume d = 2 and prove Theorem 11.11 for this case. Unlike the case of d > 3, 
due to the recurrence of planar Brownian motions, the methods in [U [20] are not applicable in 
dimension d = 2 even though we have the Harnack and boundary Harnack principles. We use a 
capacitary approach and some recent results on the subordinate killed Brownian motions instead. 

First we derive the lower bound. The method we use relies on comparing the process X a ' D , 
which is the killed subordinate Brownian motion, with another process, the subordinate killed 
Brownian motion. This method also works for dimensions d > 3. 

To be more precise, let D be a bounded open set in M. 2 and X 0,D the killed Brownian motion 
in D. Let (T t a : t > 0) be a subordinator independent of X° which can be written as T t a = t + a 2 T t 
where (Tt : t > 0) is an a/2-stable subordinator. The process (Z?' D : t > 0) defined by Z^ D = 
is called a subordinate killed Brownian motion in D. Let u a be the potential density of T a (see 
(I2.1|) ). It follows from [33] that the Green function Rf :) (x,y) of Z a,D is given by 

/>oo 

R a D (x,y) = / p° D (t,x,y)u a (t)dt, (4.1) 
J o 

where p° D {t,x,y) is the transition density of the killed Brownian motion X 0,D . It is well known 
(see, for instance, [35J Proposition 3.1]) that 

R a D (x, V) < G a D (x, V), (x, y)eDxD. (4.2) 
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Theorem 4.1 Suppose that M > 0. For any bounded C > connected open set D in M? , there 
exists a positive constant Cg = Cg(a, M, D) such that for all x,y £ D and all a S (0, M], 

G D {x,y) > R D {x,y) > C 9 log 1 + _ 

V F 2/1 

Proof. First recall the following lower bound for the transition density of the killed Brownian 
motion X 0,D obtained in [36j which states that for any A > 0, there exist positive constants cq and 
ci such that for any t £ (0, A] and any x, y € D, 

p u D (t,x,y) > c I 1 A Jt expl . (4.3) 

It follows from ((27T]) that 

u a (t) = u l (a^t) for t > 0. (4.4) 

Let T = (diam(L>)) 2 . Since u (t) is a completely monotone function with n 1 (0+) = 1, by (|4.4p . for 
any a G (0, M] 

n a (t) > w^M^T), t e (0,71. (4.5) 
By a change of variables s = — p , we have 

' 'l A WfoM) t -l e -c^ rf4 - JT^ (l A fo W fo y <f ) (46) 
Define 

Mx, y) = M^. (4.7) 
Since 1/ fD(x,y) > \x — y| 2 /diam(D) 2 = \x — y\ 2 /T, we split the last integral into two parts: 

- \ A fowfoy a s -, e - c „ & a /■ /- e - n , d \ / 1A M*wy (4 . 8) 



£ - yr / vi / V f - 2/ 



and 

■i 



s 1 (1 A(f D (x,y)s))ds > l s 1 l{ s >i/f D (x,y)}ds 

= log(/ D (x,y) VI). (4.£ 



Combining (|4T3j) and (|4T6]) - (|4T9]) . we have 

p° D (t,x,y)dt > c 2 (lAf D (x,y))+c 2 log(f D (x,y) VI) 

> C 3 logfl + ^ ( ^^ 



So it follows from (|44l)- (f42l) and (j45]) that 



G a D {x,y) > R a D (x,y) > u\M^T) j\° D (t, x,y)dt > c 4 log + 



$p{x)5p{y) 

rc - t/ M 



□ 



Integrating the estimate in Theorem 14.11 with respect to y yields the following corollary. 
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Corollary 4.2 Suppose that M > 0. For any bounded connected C ' open set D in M 2 , there 
exists a positive constant C\q = C\o(a, M, D) such that for all x G D and all a G (0, M], 

®x[td\ >C 10 5 d (x). 

The inequalities in the next lemma can be proved by elementary calculus and will be used 
several times without being mentioned explicitly. 

Lemma 4.3 For any L > 0, there exists a constant C\\ = C\\{L) > 1 such that 

C^b < log(l + b) < b for any < b < L 

and 

Cf/ log(l + s) < log(l + Ls) < C'u log(l + s) for any < s < oo. 

Using Corollary 14.21 Theorem 14.11 can be extended to general (not necessarily connected) 
bounded C ,l open sets. Recall that gf-, is defined by (|1.1[) . 

Theorem 4.4 Suppose that D is a bounded C > open set in R 2 with characteristics (R,A). There 
exists a positive constant C\i = Ci2(a, M, D) such that for all x,y G D and all a G (0, M], 

G a D (x,y)>C 12 g a D (x,y). 

Proof. Recall that y) is defined in (|4.7p . If x and y are in the same component, say x,y G U, 

then by monotonicity, 

GU*,y) > G a u(x,y) > c 1 log(l + f u (x,y)) = c 1 log(l + f D (x,y)). (4.10) 

If x,y are in the different components of D, using Corollary 14.21 and Lemma 14.31 and by following 
the second part of the proof of Theorem 11.11 in case d > 3 (that is, the paragraph containing 
(I3TT21 - (I3T3|) \ we get 

G a D (x,y) > c 2 a a 6 D {x)5 D (y) > c 2 R a a — ^— > c 3 a a log 1 H ■ ^— . 

\x — y\ z \ \x — yy J 

This completes the proof of the theorem. □ 

Recall that when d > 2 and a > 0, any non-empty open set D C M. d is Greenian for X a . For 
any Greenian open set D, any Borel subset A of D and a > 0, we define 

Cap£)(j4) := sup |?](A) : 7] is a measure supported on A 

with ^G a D (x,y) V (dy)<l}. (4.11) 

It is known (cf. [19]) that for any open subset A of D, 

C&p a D (A) = inf {£ a (u, u): ue W 1,2 (R d ),u = on D c , u > 1 a.e. on ^} 
and for any Borel subset A of D, 

Cap%(A) = inf {Cap % (B) : A C B and B is open}. 
Since £° < £ a , for any Greenian open set D C M. d and every a G [0, M] 

Cap° D {A) < Cap a D {A) for every Ac D. (4.12) 
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Lemma 4.5 There exists C13 > such that 



C13 



Cap£ (0il) (B{0, r)) > for every r G (0, 3/4). 



Proof. Recall that (see, e.g., pH p. 178]) 



G° Bm (x,y) = ^log (l + (1 '^p'^ 



Let denote the family of all probability measures on -6(0, r). It follows from p. 159] that 



-1 



Cap° (0il) (S(0,r))= inf sup G% (0>1) /i(x) . (4.14) 

x£supp(/i) / 



(4.15) 



Let m r be the normalized Lebesgue measure on £>(0,r). By (I4.14D . 

Cap° B(0 , 1)M > — ^ — . 

ieB(0,r) 

Further, by using symmetry in the first equality, and (|4.13|) in the second line, we have 
sup G° B(01) m r (x) = G B{01) m r (0) = / G° B(0 >r> (0, y) m r (dy) 

xeB(Q,r) ' JB{0,r) 

If 1 1 , 1 r 2 / , 1\ , 1 

= — 2 / 7T lo &T!2 d y = — 2 ~T 1 + 2 lo § _ < clog - , 

for some constant c > 0. This together with (14.150 yields the desired capacity estimate. □ 

For any Borel subset V, we use o~y to denote the first hitting time of V by X a : a v = inf{t > 
: X? G V}. 

Lemma 4.6 Suppose that M > 0. There exists Cm > snc/i t/iot /or even/ a G (0,M], any 
Greenian open set D in R 2 containing B(0, 1) and any x G -6(0, |) 

Cap2,(£(0,|a;|/2)) V s(o>l/2) u ) 



Proof. Fix x G 5(0,3/4) and let r := |x|/2. Since -B(0, r) is a compact subset of D, there exists 
a capacitary measure /i" for B(0,r) with respect to X a,D such that 



Gap?,(J3(0 I r))=^(B(0,r)) 
(see, for example, [21 Section VI. 4] for details). Then by Proposition 12.21 we have 



G a D (x,y)tf(dy) > inf G%(x,y) ) i£(B(0,r)) 



B(Q,r) UeB(0,r) 



> Cl G a D {x,0)Cap a D {B(0,r)) 



> Cl G a D (x,0)C a p° D (B(0,r)) (4.16) 
17 



for some constant c\ > 0. In the last inequality above, we have used (|4. 12j) . 
On the other hand, 



I G a D {x,y)n a r {dy) = f 

JB(0,r) JE 



B(0,r) 



B(0,r) 



sup 



iU P I G D (w,y)tf(dy)) F x (a^ < r D 



< 



In the last inequality above, we have used (|4.1ip . 
Combining (|4.16p - (|4.17p we have 



(4.17) 



G%(x,0) < 



Cap°,(5(M) 



B(0,r) 



□ 



Corollary 4.7 Suppose that M > 0. There exists C\$ > such that for every a G (0, M\ and 
every x G B(0, 3/4) 

G a mi) (x,0) < Cib log(l/|z|). 
Proof. It follows from Lemmas I4.5H4.6I that 



< C 1 4Cr 3 1 log(2/|i|)< C log(l/|rc|) 



for some constant c > 0. 



□ 



Lemma 4.8 Suppose that M > and that D is a bounded C > open set in M? with characteristics 
(R,A). There exists C±q = C\q(D) > such that for every a G (0, M] and all x,y G D with 
\x-y\< %6 D {x) < \R, 



G%(x,y) < Ci 6 log 



\x - y\ 



Proof. By our assumption, D satisfies the uniform exterior ball condition with radius R > 0. 

Fix x,y G D with \x — y\ < |5d(x) and let r := 5d(x). Since r < R, without loss of generality, 
we may assume x = (0, 0), (0, -r) G dD with B((0, -2r),r) G R 2 \ D. 

Let a := ar (2 - a V a , y := r~ x y and D := r~ l D. Then by <|212|> . 



G^(0,y)=G|(0,y). 
By the strong Markov property, we have 



G|(0,2/)=G| (0il) (0,y)+E 



U S(0,1) 



(4.18) 



(4.19) 
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Note that a = a A 2 - a ^ a < MR {2 ~ a "i/ a . Define 

h(z, w) := E z 



G%(X% ,w) 

U S(0,1) 



For each fixed z G B(0, 1), the function w \— > h(z, w) is harmonic in B(0, 1) with respect to X a and 
for each fixed w G B(0, 1), z t— > h(z,w) is harmonic in 5(0, 1) with respect to process X a . So it 
follows from Proposition 12,21 

h(0,y) < ci min h(z,w) < c\ min G fl -(z,«i) < ciG~(0, x±) 
where |xi| = 1/2. In the second inequality we used that G^(-,w) is superharmonic in B(0, 1) for 

P 2 



X a . Note that Dc£:=l \ B((0, -2), 1). Thus by Lemma 



_C2_ 

Cap^(B(0,l/4)) 



/i(0, £) < Cl G a s (0, Xl ) < ciG|(0, Xl ) < - n/ J - < oo. (4.20) 



On the other hand, by Corollary 14.71 

G| (0)1) (0, y) < c 3 log = c 3 log . (4.21) 

It follows from KTEh - KTD i that 

GUx, y) = GUO, V)<c 4 + c 3 log (^) < c 5 log 

V \y\ J \\ x -y\ 

which proves the lemma. □ 

Lemma 4.9 Suppose that M > and that D is a bounded C 1, open set in M?. If x and y are in 
the same component of D with 

-(6 D (x) V S D (y)) <\x-y\< c{5 D {x) A S D (y)) 
c 

for some c > 1, then there exists Cyj = Cu(c, D) > such that for every a G [0, M] 

G a D {x,y) < Ci 7 . 

Proof. Without loss of generality, we assume Sd(x) < Su(y). If \5oiy) < |# — y| < ^5d(%), then 
the lemma follows from Lemma 14.81 In the case jSd(x) < \x — y\ < c5d(x), since x,y are in the 
same component of D, we use Proposition 12.21 and a standard Harnack chain argument. □ 

Theorem 4.10 Suppose that M > and that D is a bounded C 1 ' 1 open set in M 2 . There exists 
C±s = Cis(D) > such that for every a G (0, M\ and all x,y G D 

G a D (x,y) < C 18 log(l + 



x — y\ 2 
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Proof. Let L = max{2diam(D),2}. 

(i) If \x — y\ < 1/4, by Lemma 14.81 applied to B(x,L), 

Gb(?,v) < G a B{x , L) (x,y) < Cl log (|^) < c 2 log (l + 



(ii) If 1/4 <\x - y\, by (pl2|) and Corollary ST 



G%(x,y) < G% {XtL) (x, y) < G^Z^L^x, L~ l y) < c 3 log ' L 

1 



\x - y 



< c 3 log (4L) < c 4 < c 5 log 1 + 



\x — y\ 2 



□ 



Now we are ready to prove Theorem 11.11 for d = 2. Recall that fo(x,y) is defined in ()4.7[) . 

Proof of Theorem ll.il when d = 2: By Theorem 14.41 we only need to consider the upper bound. 
We divide its proof into two steps. 

Step 1. We first consider the case that x and y are in the same component of D. Without loss 
of generality, throughout this proof, we assume 5d(x) < &d{v)- 

Fix zq £ D with kR < 8d(zq) < R, and let e± := kR/2A. Choose Q x , Q y £ 3D with \Q X — x\ = 
6d(x) and \Q y — y\ = 5o(y)- We consider the following five cases separately. 

(a) If 5 D (x) > £ik 2 /32, by Theorem SITO] 

G%(x,y) < cilog^l+ - c 2 log(l + / D (x,y)). 

(b) Suppose 5d{x) < £ik 2 /32 and 5d{u) > E\k/A. Let r := £ik/16 and put x\ = A rK / 2 {Q x )- One 
can easily check that \zq — Q x \ > r and \y — Q x \ > r. So by (|2.9p . Theorem 12.31 and Theorem 14. 101 
we have 

G a D (x,y) < c 3 G a D ( Xl ,y) ^ X,Z0 \ < c A 5 D {x) < c 5 f D (x,y) 

for some 03,04,05 > 0. Note that fo{x,y) < cq in this case because D is bounded and \x — y\ > 
^d(v) — $d(x) > £ik(1/4 — k/32) > 0. So it follows from the above display and Lemma 14.31 that 

G a D (x,y) < c 7 log(l + f D (x,y)). 



(c) Suppose 8d(x) < £ik 2 /32, Sd(v) < e\k/A and \x — y\ < ^£)(y)/2. From 5d{v) < \x — y\ + 5d{x) 
we conclude that &d{v) < 25d(x) and so \x — y\ < 5d(x). This together with Lemma I4~8l gives that 

G%(x,y) < cslog(^L) < c 9 log(l + / D (x,y)). 

\\ x -y\J 

(d) If \8 D (y) <\x-y\< (24/k 2 )5 d (x), by LemmaSSl 

G a D {x,y) < cio < cnlog(l + f D (x,y)). 
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(e) The remaining case is 



S D (x) < 8 D (y) < ^ and |ac - y| > max j ^<5 D (x), 



We claim that in this case 

G%(x,y) < c 12 f D (x,y). (4.22) 
By Lemma I4.3| the above implies that G a D (x, y) < C13 log(l + /d(sc, y)) since in this case 

f ~~ v\ 

We now proceed to prove (|4.22p by considering the following two subcases. 

(i) (24/ 'k 2 )5d(x) < \x — y\ < {&/ k)5d{v)'- Let r := 5o(y)/3. Put x\ = A rK / 2 (Q x )- One can easily 
check that \zq — Q x \ > r and \y — Q x \ > r. So by (|2.9|) and Theorem 12.31 we have 



Moreover, 



implying that 



24 3 6 72 

—5 D (x) <\x-y\< -\xi - y| < -S D (y) < -^(rci), 



36 N 

Sd(v) < \x - y\ + 5 D (x) < I ^ + -x ) |xi - y\ . 



It follows from Lemma 14.91 that Gp(xi,y) < c±q. Therefore 

G a D {x,y) < C17 f D X , < cisfD(x,y). 
\x - y\ 

(ii) <5d(z) < <5 D (y) < (w/4)|a;-y|: Let r := ^(\x-y\ Aei). Put zi = Ave/2 (Qa:) and yi = A rK / 2 (Qy). 
Then, since |zo — Qzl > r and |y — Q x | > r, by (|2.9p . we have 

_! G£(si,y) G^(x,y) . G^(x l5 y) 

c 19 T^TT- TT < 7^ 77 — ^ C19: 



GS,(xi,^o) - G^(rr,z ) " 13 G% (zi, 2 ) 
for some C19 > 1. On the other hand, since \zo — Q y \ > r and |xi — Q y | > r, applying (|2.9p . 

_! G a D {xuVi) < G a D (yi,*o) < G a D (x uyi ) 
Cl9 G^y) " G-(y,z ) " ° 19 G^, y) " 

Putting the four inequalities above together we get 

G a D ( Xl , yi ) < G^(rr,y) < c2 G a D ( Xl , yi ) 



19 Gf^^GMyi^o) " ^(x,^ )G^(y,z ) ~ 19 G£>K zo)G a D ( yi , z ) ' 
Moreover, ||x — y| < |^i — yi| < 2|x — y| and 

4 1 64 

— (fe(xi) V 5 D (yi)) < -\x - y\ < \x\ - y x \ < -3— OM^i) A fo(yi))- 
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Thus by Lemma 14.91 and Theorem 12.31 we have 



m f \^ G a D (x,z )G a D (y,z ) 

G D {x,y) < C20 — - <c 2 if D {x,y) 

G D {Xi, z 0)G D {yi, z 0) 

for some C2o,C2i > 0. This completes the proof of the claim (|4.22p and therefore of the theorem 
when x and y are are in the same component of D. 

Step 2. Next we consider the case that x and y are in two different components of D. This 
part of the proof is the same as the second part of the proof of Theorem 11.11 when d > 3 (that 
is, the paragraph containing (|3.12p - (|3.13p ). The only place that needs modification is the proof of 
sup zeD E x [t£)] < C22 < oo. When d = 2, we can not use ()3.5[) to deduce it. However, since D is 
bounded, there is K > so that D C B(0, K). It follows from Step 1 that 

sup E^[r^] < sup ^z[t% { q K) ] < c 23 < oo. 

zeD,ae(0,M] zeB(0,K),ae(0,M] 

This completes the proof of Theorem II. H □ 



Theorem 4.11 (3G theorem for d = 2) Suppose that M > and that D is a bounded C 1, open 
set in R 2 . Then there exist positive constants C\g = C\g(D, a, M) and C20 = C2o{D,a,M) such 
that for all x,y, z G D and a € (0, M] 

GU^y)GUy,z) < c 19 (io g (i + / D ( 2; , y )) + iog(i + / D ( y ,.)) + i) 

Lr D {X, Z) 

< C 20 f (log^^ VI) + (log^^ Vl) j . 

\ \x — y\ ' \ \y — z\ J J 

Proof. Note that, if x,z are in different components of D, either x,y or y,z are in different 
components of D. Thus, by Theorem II. II for d = 2 and the fact that a G [0, M], we have 

G a D (x,y)G a D (y,z) log(l + f D (x, y)) log(l + f D (y, z)) 
G a D (x,z) ~ Cl log(l + f D (x,z)) 

for some c\ = c\(M,D,a). Now following the proof of |16^ Theorem 6.24], we get the theorem. □ 



Remark 4.12 By considering how many different components of D that x, y and z fall into, we 
could get more precise 3G estimates with the dependence on a explicitly spelled out. Theorem l4.11l 
will not be used in the remainder of this paper. 

5 One Dimensional Case 

In this section we assume d = 1 and prove Theorem 11.11 for this case. We will follow the ideas in 

Let X be the supremum process of X a defined by X t = sup{0 V X s : < s < t} and let 
X~ a — X a be the reflected process at the supremum. The local time at zero of X a — X a is denoted 
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by L a = (Lf : t > 0) and the inverse local time by {rf : t > 0}, where r° := inf{s : > t}. 
The inverse local time {r° : t > 0} is a subordinator. The (ascending) ladder height process of 
X a is the process H a = (Hf : t > 0) defined by = X® a . The ladder height process is again a 
subordinator. It follows from [26] that iif a is a special subordinator with Laplace exponent given 
by 

X (A) = exp ^- ^ dflj (5.1) 

and that the drift coefficient of H a is 1. When a = 0, we have x°(A) = A. Thus, if V a is the 
potential measure of H a and V a (x) = U a ([0, x]), then, for every a > 0, V™ has a continuous, 
decreasing and strictly positive potential density v a such that v a (0+) = 1. When a = 0, we have 
v a = 1. The following results is a uniform version of [261 Proposition 2.3] in our present special 
case. 

Lemma 5.1 Let M and R2 be positive constants. There exists a constant C21 = C2\{M, R2) G 
(0, 1) such that for all a G [0, M] and x G (0, R2], 

C21 < v a (x) < C^i and C 2 ix < V a (x) < C^x . 

Proof. Since H a is special, the potential density v a is a decreasing function. Hence info<t<_R 2 v a (t) = 
v a (R2). It follows from (15.11) that the Laplace exponent x a is continuous in a. Thus, the potential 
measures converge vaguely, and by continuity and monotonicity of v a , we get that v a (t) — ► v b (t) 
as a —>■ b for all t > 0. In particular, v a (R2) — > v b (R2)- Therefore c\ := info<t<_R 2 ,o<a<M v a (t) > 0. 
Since v a (t) < 1 for all t > and all a > 0, we get that C2 = sup 0<4<J j 2 o< a <M v a (t) = 1. Choose 
C3 = C3(M, R2) G (0,1) such that C3 < ci < C2 < cj 1 . Since V(a:) = fyv(t)dt, the claim follows 
immediately. □ 



Theorem 5.2 Suppose that M > 0. For any bounded open interval D in R, i/iere exists a constant 
C22 = C22(a, M, D) > 1 suc/i i/iat /or all x,y G D and aZZ a G (0, M], 

n -\ f rx i \ x ( \a1/2 a M^0My)\ , na , \ ^ r 1 ( fx 1 \x 1 ^^l/2 A MgMgtj/) 
C 22 ((5i)(x)(5 D (y)) 7 A — : : — < G D {x,y) < C22 {8 D (x)5 D {y)) 1 A — ; — 



x — y\ J \ \x — y 

Proof. The proof of the lower bound is similar to that of Theorem 14.11 and [261 Proposition 3.3]. 
Using our Lemma 15.11 instead of [26|. Proposition 2.3], we can follow the proof of [261 Propsoition 
3.1] to get the upper bound. We omit the details. □ 

Integrating the estimate in Theorem 15.21 with respect to y yields the following corollary. 

Corollary 5.3 Suppose that M > 0. For any bounded open interval D in R, there exists a positive 
constant C23 = 623(0:, M, D) > 1 such that for all x G D and all a G (0, M), 

C^S D (x) < EM] < C 2 ^ D {x) . 

Using Corollary 15.31 we can repeat the proof of Theorem 11.11 for d > 3 case (see also [26\ 
Theorem 3.8]) to generalize Theorem 15.21 to general (not necessarily connected) bounded C 1 ' open 
sets. Recall that g a D is defined by (|1.2|) . 
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Theorem 5.4 Suppose that D is a bounded C > open set in K. There exists a positive constant 
C24 = 624(0;, M, D) > 1 suc/i £/iat /or aZZ x,y £ D and all a G (0, M], 

C^g D {x,y) < G a D (x,y) < C 2A g a D {x,y). 



6 Martin Boundary and Martin Kernel Estimates 

Throughout this section we assume that d > 1 and D is a bounded C 1 ' 1 open set in R rf with 
characteristics (R, A) and the corresponding k. We will show in this section that the Martin 
boundary and the minimal Martin boundary of D with respect to X a can both be identified with 
the Euclidean boundary dD of D. With the boundary Harnack principle given in Theorem l2.3l the 
arguments of this section are modifications of the corresponding parts of [31 [T4"l [23l [25] . For this 
reason, most of the proofs in this section will be omitted. 
The next lemma follows from Theorem 12.31 

Lemma 6.1 Suppose that M > and that D is a bounded C 1 ' 1 open set in M. d . There exists a 
positive constant C25 = C25(D,a, M) such that for all a G (0, M], Q G dD, r G (0,R/2), and 
nonnegative function u in Mr which is harmonic with respect to X a in D n B(Q,r) we have 

k 



u{MQ)) < C25 ( - ) u(A {K/2)kr (Q)), k = 0,l 



(6.1) 



Lemma 6.2 Suppose that M > 0. For every b G (0,oo), there exist C2& = C2e(M,b) > and 
C27 = C 27 (M, b) > such that for all x Q G R d , a G (0, M] and r G (0, 6], 



C 26 r 2 < E, 



and 



r S(x ,r) 



< C 27 r 2 



(6.2) 
(6.3) 



Proof. See |32[ Lemmas 2.2 and 2.3] or |13[ Lemmas 2.3 and 2.4] for a proof of (|6.2p . The inequality 
()6.3p follows easily from Theorem ll.il In fact, by (|2.1ip and Theorem 1 1 . 1 1 (with Mb^ 2 ~ a ^ a instead 
of M) 



E., 



B(0,r) 



B(0,1) 



r B(0 



1} (r x,z)(iz < cr ^(0,1) a; ) = cr 5 B{0 , r) (x). 



□ 



For an open set U C M d , let 



tfJKx.z) := / G^(x,y)J a (y,z)dy, {x,z)eUxW. 
Ju 

Then by (12. 4p . for any non- negative measurable function / on M. d , 



(6.4) 



E, 



f{ X T? T )'l X T? T - / ^T? r 



^(x,z)/(z)dz. 



From (|6.4p . Theorem 1 1 . 1 1 and Lemma 16.21 we immediately get the following proposition. 



24 



Proposition 6.3 Suppose that M > 0. There exist C 2 s > and C29 > such that for all 
a £ (0, M] and r £ (0, R) and x £ R d , 



K a BM (x,y) < C 2 8r(r-\x-x \)(\y-xo\-r)- d - a for (x,y) £ B(x ,r) x B(x ,r) (6.5) 
and 

K B( xo ,r)(xo,y) > C 29 r 2 \y-x \- d - a foryeB{x , r y. (6.6) 

Using (|6.6p . the proof of the next lemma is similar to that of |23^ Lemma 4.3] or [25|, Lemma 
5.3]. Thus we skip the proof. 

Lemma 6.4 Suppose that a > and that D is a bounded C 1,1 open set in M. d . There exists a 
positive constant C30 = Cso(D, a, a) such that for all Q £ dD, r £ (0, R/2) and w £ D\ B(Q, r), 

G a D (A r (Q),w) > C 30 r 2 [ f(\z-Q\/2)G a D (z,w)dz. 

JB{Q,rY 

Using (|6.5p , Lemmas 16.11 and 16.41 the proof of the next lemma is similar to that of |23} Lemma 
4.4] or |25l Lemma 5.4]. Thus we skip the proof. 

Lemma 6.5 Suppose that a > and that D is a bounded C 1,1 open set in M. d . There exist positive 
constants C31 = Cs±(D,a, a) and C32 = C32(D,a,a) < 1 such that for any Q £ dD, r £ (0, R/4) 
and w £ D \ B(Q, 2r/n), we have 



G a D (X% B ,w):X% £B(Q,r) c < C 31 C^G a D (x,w), x£DnB k , 



where B k := B(Q, (n/2) k r), k = 0, 1, ... . 
Let xq £ D be fixed and set 

M a D (x,y) := ^ X ' y) x,y£D, y + x . 

Now the next theorem follows from Theorem 12.31 and Lemma 16.51 (instead of [3l Lemma 13] 
and [31 Lemma 14], respectively) in very much the same way as in the case of symmetric stable 
processes in [3l Lemma 16] (with Green functions instead of harmonic functions). We omit the 
details. 

Theorem 6.6 Suppose that a > and that D is a bounded C 1,1 open set in M. d . There exist positive 
constants R\, M\, C33 and (5 depending on D, a and a such that for any Q £ dD, r < R\ and 
z £ D\ B(Q, M x r), we have 

\MUz,x)-MUz,y)\ < C33 (t^y[\ x,y£DnB(Q,r). 



In particular, the limit lim.£) B y^ w Mf ) (x, y) exists for every w £ dD. 
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As the process X a,D satisfies Hypothesis (B) in Kunita and Watanabe [27], the process X a,D 
has a Martin boundary: For every a G (0, M], there is a compactification D^f of D, unique up 
to a homeomorphism, such that Mf)(x,y) has a continuous extension to D x (D*f \ {xq}) and 
Afg(-,zi) = M£(-,z 2 ) if and only if zi = z 2 (see, for instance, [27]). The set d^D = Djf\D is 
called the Martin boundary of D for X a ' D . For z G c^A set M£(-, z) to be zero in Z) c . 

For each fixed z G dD and x G -D, let 

MSfoz) := lim Mf,(x,y), 

DBy^z 

which exists by Theorem 16.61 Mp is called the Martin kernel of D with respect to X a . For each 
z G dD, set Mf)(x,z) to be zero for x G D c . By Theorem 16.61 Mp(z,x) is jointly continuous on 
{x e D : 5d(x) > 2e} x dD, and hence on D x 51? after letting s J. 0. 

The following Martin kernel estimate is an immediate consequence of the Green function esti- 
mates in Theorem I l.li Recall that D{x) denotes the component of D that contains x. Define 

if x G D(x Q ),z G 3D(x ) or x G D \ D(x ), z£dD\ (dD(x ) U &D(x)), 
if x G £>(x ), z£dD\ dD(x ), 
if x G D\D(x ),z G <9L>(x ). 

(6.7) 

Theorem 6.7 Suppose that M > anci t/iai D is a bounded C 1 ' 1 open sei in There exists 
C34 := C34(xo, -D, a, M) > 1 suc/i i/iaf /or a// a G (0, M], 

Of/ h a D {x, z) < M%{x, z) < C 34 h a D (x, z) forxeD,ze dD. 



h a D (x,z) :-- 



( Sp{x) 



\x-z\ d 
a S D (x) 
\x— z\ d 



Theorem 16.71 in particular implies that Mfy(-,z\) differs from Mp^z?) if Z\ and 22 are two 
different points on dD. 

Now using our Green function estimates, (|6.5p and Lemma 16.14 one can follow the arguments 
in the proofs of [23\ Lemmas 4.6-4.7] and \25\ Lemmas 5.6-5.7] and get the next two lemmas. 

Lemma 6.8 Suppose that D is a bounded C 1 ' 1 open set in M. d . For every z G dD and B C B C D, 
M%{X"a,z) isP x -integrable. 

Lemma 6.9 Suppose that D is a bounded C 1 ' 1 open set in M, d . For every z G dD and x G D, 

M%(x,z) = E x \m%(X^L ,z)] for every V <r <\{Rf\5 D {x)). (6.8) 

Unlike the case in the proofs of [3 Theorem 2.2] and [23j Theorem 4.8], P x {X? a G dU) / 
for every smooth open set U . Thus we give the details of the proof of the next result. 

Theorem 6.10 For every z G dD, the function x \— > M^(x,z) is harmonic in D with respect to 
X a . 
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Proof. Fix z G dD and let h(x) := M^(x,z). Consider an open set D± C D\ C D and x G -t>i 
and put 

r(x) = -(R A 5d(x)) and B(x) = B(x,r(x)). 
Define a sequence of stopping times {T m ,m > 1} as follows: 

Ti = inf{t >0:X?$ B(Xg)}, 

and for m > 2, 

Tm.-1 + Tspei ) O 6»T m _i 
v T m-1 



L m- 
r" 



if X£ G Di 

^ m— 1 

otherwise. 



Note that X"^ G dD\ on n^^T^ < rj^}. Thus, since lrm m ^. oo r m = P x -a.s. and /i is 
continuous in l), 

lim fr(X£J = h{X a T a ), on n- 1 {T n < rgj 
and, since /i is bounded on by the dominated convergence theorem 

lim E x [h(X£J; n» {T n < r^}] = E x ); n^ = {T n < r^} . 

Therefore, using Lemma 16.91 
h{x) = lim E* [h(X£J] 

m— >oo L J 

= lim E x [/*(X£J; U~ {T n = T a D J] + lim [fc(*£J; fl~ {T n < rgj] 



E, 



E, 



); U~ {T„ = t£J +E, ); n^ = {T n < rgj 



□ 



A consequence of Theorems 16. 6| 16.71 and I6.1UI is that, when D is a bounded C > open set, the 
Martin boundary of X a ' D can be identified with the Euclidean boundary dD of D. 

A positive harmonic function u for X a,D is minimal if, whenever v is a positive harmonic function 
for X a ' D with v < u on D, one must have u = cv for some constant c. The set of points z G d^D 
such that M^(-,z) is minimal harmonic for X a,D is called the minimal Martin boundary of D for 
X a > D . 

With the explicit estimates from Theorem 16.71 by the same argument as that for |14l Theorem 
3.7], we have the following. 

Theorem 6.11 Suppose that D is a bounded C 1,1 open set in W 1 and a > 0. For every z G 3D, 
Mp(-,z) is a minimal harmonic function for X a,D . Thus the minimal Martin boundary of D can 
be identified with the Euclidean boundary. 

We know from the general theory in Kunita and Watanabe [27] that non-negative superharmonic 
functions with respect to X a,D (or equivalently, superharmonic functions with respect to X a that 
vanish on D°) admit a Martin representation. Thus, by Theorem 16. 1 II we conclude that, for every 
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super harmonic function u > with respect to X a ' , there is a unique Radon measure /i in D and 
a finite measure v on <9D such that 

u(x) = [ G a D (*,y)Kdy)+ [ M%(x,z)v{dz). (6.9) 

JD JdD 
Furthermore, u is harmonic for X a,D if and only if the measure fi = 0. 

7 Perturbation Results 

In this section we assume d > 1 and fix a > 0. We consider a symmetric Levy process Z which can 
be thought of as a perturbation of X a , and show that under certain conditions, the Green function 
of Z D , the process Z killed upon exiting a bounded C 1,1 open set D, is comparable to the Green 
function of X a,D , see Theorem 17.131 Together with Theorem [TTTJ this gives sharp bounds for the 
Green function of Z D . 

The approach of this section is motivated by [T7], where perturbations of pure jump Levy 
processes are discussed. Even though they consider pure jump Levy processes, some results work 
for our case as well. 

Throughout this section, Z is a symmetric Levy process in M. d such that its Levy measure has 
a density J z (x, y) = j z (y — x) with respect to the Lebesgue measure. We assume that 

#(s) :=f(x)-j z (x) 

is nonnegative and integrable in W d and put J a := J^ d jf(y)dy. We also assume that the transition 
density of the Levy process Z exists and we denote it by p z (t, x, y) = p z (t, y — x). 

Recall that p a (t,x,y) = p a (t,y — x) is the transition density function of X a . It is well known 
that (see (HI El]) 

p a (t, x,y) < c (V d/Q A t~ d/2 ) A (V<V2 e -c2|*- y | 2 A + * ^ , (t, x, y) £ (0, oo)xR d x R d 

(7.1) 

for some c = c(a, d, a) > 0. Thus, by following the proof of [TTJ Lemma 2.6], we have the following. 

Lemma 7.1 p z (t,x) is bounded on the set {(t,x) :t>0,\x\ > e} for e > 0. 

Recall that for every bounded open set D, G a D is the Green function of X a ' D . We know from 
[244 Corollary 3.12] that there is a constant c = c(D,a) such that 

cE x [t%] E y [t%] < G a D (x, y), x,yeD. (7.2) 

The proofs of the following three results are the same as those of [171 Lemmas 2.2, 2.4, 2.5]. So 
we omit their proofs here. In the remainder of this section, the dependence of the constants on Z 
will not be mentioned explicitly. 

Lemma 7.2 For every bounded open set D, there exists C35 = C^{D,a) > such that for all 
x G D and t > 1 we have 

p a D (t,x,y) < C^r 2 E x [T a D ] E y [T a D ]. 
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For any open set U C M. d , r z := inf{i > : Zt ^ U} denotes the first exit time from U by Z. 
We denote by Z D the subprocess of Z killed upon leaving D and pf)(t,x,y) the transition density 
for Z D . 

Lemma 7.3 For every bounded open set D, there exist a constant C^q = Cz§{D,a) > 1 such that 
for every x £ D, 

C^MtE] < E x [rh] < C 3G E x [rE). 
Lemma 7.4 For every bounded open set D and any x £ D and t > 0, we have 

pUt,x,-)<e Jat p%(t,x,-) a.e. 
If, in addition, we assume that p z (t, •) is continuous then we have for x,y £ D, 



E, 



p z (t-TE,Z r z,y) :*>rg 



< e 



2 J at 



p a (t-T a D ,X%,y):t>r a D 



(7.3) 



Using the above lemmas (for Lemma 17.41 only its first part is needed), and following the proof 
of [171 Theorem 3.1], we have 

Theorem 7.5 For every bounded open set D, there exists C37 = C^{D^a) > such that for every 
x £ D, 

G z D (x,y) < C 37 G a D (x,y) a.e. y £ D. (7.4) 

Here are some assumptions that we might put on the process Z. 

(Al) The transition density pfj(t, x, y) of Z D is continuous and strictly positive in (0, 00) x D x D. 
(A2) There exist positive constants c and p such that jf(x) < c\x\ p ~ d on B(0, 1). 

(A3) j z satisfies either (a) or (b) below: 

(a) There exists a non-negative Borel function L(x), which is locally integrable on {0}, 
such that for any Borel set B, \B\ = J B L(x)j z (x)dx. 

(b) There exists Rq > such that inf^g^o^o) j z {x) > 0. 

Without loss of generality, we may and do assume that the constant p is less than 1. 

Proposition 7.6 ([24], Corollary 3.11]) Suppose that (Al) and (A3) hold. Then for every bounded 
open set D, there exists constant Cs$ = Css(D,a) > such that 



C 38 E x [T Z ]E y [T Z ] < G z D {x,y), for all (x,y) £ D x D. 
For the remainder part of this section, we assume D is a bounded C 1,1 open set in M. d . 



(7.5) 



Lemma 7.7 Suppose that (Al) and (A3) hold. Then for every e > 0, there exists C 39 = C3g(e, D,a) > 
such that for all x,y £ D satisfying \x — y\ > e, 

G a D (x,y) < C 39 G z D (x,y). 
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Proof. It follows from Theorem 11.11 that there exists c\ = c\(D,a) > such that 

cT/X [rg] < 5 D (x) < Cl E x [T a D ] , x G D. (7.6) 

Combining (|7.6p with Theorem 11.11 yields that there exist C2 , C3 > so that for all x and y G D 
with |x — y| > e, 

G a D (x,v) < c 2 fe(x),5 D (y) < c z E x [r%}E y [r%}. 
Therefore, by Lemma 17.31 and Proposition 17.61 we get 

G a D (x,y) < cA[rE}E y [TE] < c 5 G D (x,y) 

for some positive constants 04,05 > 0. □ 

The next lemma can be proved by following the arguments in the proofs of [311 Lemmas 7, 9]. 
So we skip the details here. 

Lemma 7.8 Suppose that (Al) holds. For all x,w G D, we have 



G a D (x,w) < Gf,(x,w) + / / G a D (x,y)j?(y-z)G a D (z,w)dydz. 

J D J D 

Proposition 7.9 Suppose that d > 3 and that (A2) holds. There exists a positive constant C40 
Cao{D,o) such that for every (x,w) G D x D 




G a D (x,y)j a (y - z)G a D (z,w)dydz < C m G a D (x, w)\x - w\ d ~ 



Proof. Using the generalized 3G inequality (Theorems I3.5|) for the Green function of X a , one can 
easily get the following 

Gh(?,y)ft(y - z)G a D {^w)dydz 

V — z\ p ~ d dydz ^ d-1 II \y — z\ p ~ d dydz 



'D JD 

< cG D (x,w)(\x-M d - 2 [ [ | 'V 2 , 2 + \x-w\ d ^ 

V JDJD\x-y\ Z \z-W\ d 1 JDJD\X~yr X \Z-W 

\y — z\ p ~ d dydz id f f \y ~ z\ p ~ d dydz 



+ \x-w\ d 1 / / | — - , , ' , ".,-. + \x-w, 

I D J D \ x ~ V\ \Z — W\ a 1 JDJD\ x ~y\ \Z — W\ a 1 

for some constant c = c(D,a) > 0. Now combining the above with [17^ Lemma 3. 12], we easily get 
the conclusion of the proposition. □ 



Proposition 7.10 Suppose that d > 1 and that (A2) holds. There exists a positive constant 
C41 = Gn(D, a) such that for every (x,w) G D x D 




G D {x,y)j a (y - z)G D {z,w)dydz < C, 



8d(x)5 d (w) 



41 I \d-o 

IdJd \x-w\ d p 
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Proof. Recall that for every d > 1, G a D (x,y) < ci^j^# by Theorem O Therefore, by 
following the arguments in the proof of \31\ Lemma 8] we have 

/ G D (x,y)- *dy<C2 - (7.7) 

Thus 

8 D (x) 



G a D (x,y)j?(y-z)d y <c 3] 



In ~ ' " \x - z\ d ~P 

and so, since G D (z,w) = G D (w,z), by (17.7j) 



/ / G D (x,y)j a l {y-z)G D (z,w)dydz < c 3 5 D (x) [ G a D (z,w)- \j-dz 

JD JD JD \ x ~ z \ H 



< C4- 



\x — w\ d P 

□ 



Lemma 7.11 Suppose that d G {1,2}, T > and i/iai (Al) holds. Then there exists a constant 
C42 = C*42(a, T) > s-uc/t £/iai 

sup (p a (t,x,y)-e- 2J « t p Z (t,x,y)) < C^i 1 "*/ 2 . 
0<t<T 

Proof. Recall that p z (t, y — x) := p z (t, x, y) and p a (t, y — x) = p a (t, x, y). Let Z be a pure jump 
Levy process with Levy density jf in M. d independent of Z. Then Z is a compound Poisson process 
with transition probability given by 



n=l 



The process Z + Z has the same distribution as X a . Thus the distribution of X£ is equal to the 
convolution of p (£, •) and P z (t, •). Consequently, we have 



00 +n^Z 



*V(V)* (if )*"(*) 



n=l 



It follows from Lemma 17.41 and (|7.ip that for < t < T 

P Z (t, •) * Ol )*"(*) < e^V (*, •) * (tf < c ie Jat (J a ) n (t- d ^ A ^ 2 ) < c 2 t- d ^ t (J- a ) n 

for some positive constants ci, C2. Thus it follows from Lemma [7.41 and the above two displays that 
for < t < T 

p a (t, x) - e~ 2Jat p z (t, x) = p a (t, x) - e~ Jat p z (t, x) + e~ Ja \l - e~ Jat )p z (t, x) 

00 j.n—d/2( q \n 

< c 2 y- ^L + (l_e-^V(^)- 

n=l 

bmce by (HU) p a (t, x) < c 3 t~ d / 2 for < t < T, we reach the conclusion of the lemma in view of the 
above display. □ 

We also need the following lemma. 
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Lemma 7.12 Let D be a C > open set with C ' characteristics (R,A). Then there is a constant 
C43 > such that for all x G D and t > 0, 

r x (r a D >t) < C 43 ^4- 
t + v t 

Proof. When t > 1, the above inequality follows immediately from Markov's inequality and (|7.6p . 
To establish the inequality for the case of < t < 1, we will use a result from |lUj . 

We will only give the proof for the case d > 2. The proof in the case d = 1 is similar but 
simpler. Without loss of generality, we can always assume that R < 1 and A > 1. By definition, for 
every Q G dD, there is a C^-function cf> Q : R d ~ l R satisfying Q (O) = 0, V0q(O) = (0, . . . , 0), 
II V0q || oo < A, \V(j)Q(x) — V(J)q(z)\ < A\x — z\, and an orthonormal coordinate system CSq: 
y = (y, y d ) such that B(Q, R)nU = {y = (y, y d ) G B(0, R) in CSq : y d > <f>{y)}. Define 

p Q (x) := x d - 4>q{x), 

where (x, x d ) are the coordinates of x in CSq. Note that for every Q G dU and x G B(Q, R) n U, 
we have (1 + A 2 ) -1 / 2 pq(x) < 5u{x) < pq{x). We define for ri,r2 > 

DQ(r 1: r 2 ) := {y £ll :r 1 > p Q (y) > 0, \y\ < r 2 } . 

Note that for b > 0, 

P,(t£>1) < ^(^o,n0 >1 )+ P -(^ Q( , , ro) eDand ^^ 



< E x 



T D Q (5 ,r ) 



+ F X [XK ED 



T DQ(S Q ,r Q ) 

Thus by |1UI Lemma 3.5], there is a constant c\ = ci(R, A, a) so that for every b G (0, a] 

T x {T b D > 1) < cife(x) for every x G D. (7.8) 

Note that for < A < 1, X^D is a C ' open set with C > characteristics (R, A). Hence by the 
scaling property of X a in (|2.10p . we have from (|7.8p that for t G (0, 1], 



x (r a D >t) = ¥ t - 1/2 lr^ 1/2D " >1 < Cl 6 t - 1/2D (t-^x) = Cl 



This completes the proof of the lemma. □ 

Theorem 7.13 Suppose that (Al)-(A3) hold. There exists C44 = C44(D,a) > such that 

CllG z D {x,w) < G a D {x,w) < C u Gf){x,w), (x,w) G D x D. (7.9) 

Proof. By (|7.4p and Lemma [7.7l we only need to show the second inequality in (|7.9p for \x— y\ 2 < e, 
where e G (0, 1) is a constant to be chosen later. We consider the cases d > 3 and d < 2 separately. 

(a) d > 3: Applying Lemma [7T81 and then Proposition 17.91 we get 

G a D (x,y) < G z D (x,y) + Cl G^(x,y)|x-2/| d - 2 . 
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Choose e > small so that 



c 1 G a D (x,y)\x-y\< 1 - 2 < ^G a D (x,y) if \x - y\ < e 1 ' 2 . 



Thus 



G a D (x,y) < 2G z D (x,y) if \x - y\ < e 



1/2 



(b) d < 2: We first note that, since p D (t, x, y) < C2(t d l a A t d / 2 ) by (|7,ip . using the semigroup 
property, 



p a D (t,x,y) = / p a D (t/3,x,z) / p a D (t/3,z,w)p a D (t/3,w,y)dwdz 
Jd Jd 

< c 2 (t~ d/a M~ d / 2 ) [ p a D (t/3,x,z)dz [ p a D (t/3,w,y)dw 
Jd Jd 

= c 2 {t- d ' a A r d ' 2 ) F x (t% > t/3)¥ y {r a D > t/3). 

By Lemma \7.12\ we get 

Ph(t, x, y) < c 3 (t-*'° A t- d ' 2 ) Al)(^A l) 

< c 4 (i-*/*- 1 A t~ d ' 2 ~ l ) 5 D {x)5 D {y). 

Consequently, 

f 00 

/ p D (t,x,y) dt < c§t &D{x)$D(y) for every to > and x,y E D. 
Jt 

Since 

pfc(t, x, y) = p a (t, x, y) - E x [p a (t - T a D ,X^y) : t > r£ 

and 

p D (t, x, y) = p z (t, x, y) - E x [p z (t - r|, Z T z,y) : t > 
it follows from (17.31) that 



(7.10) 



p D (t,x,y) = p a (t,x,y)-E x \p a (t - X*«,y) : t > r£ 

< p a (t, x, y) - e~ 2Jat E x [p z (t - t z , Z T z,y):t> r z 

= p a {t, x, y) - e~ 2Jat (p z (t, x, y) - p z D {t, x, y)) 

< p a (t, x, y) + pf>(t, x, y) - e- 2Jat p z {t, x, y). 

So integrating over [0, to] with to = (£d(x)<5d(?/)) 1//2 , which is bounded by diam(D), we have by 
Lemma [7TE] and (TTlOl) that 

G a D (x,y) = / P a D (t,x,y)dt+ / p a D (t,x,y)dt 

JO Jt 
< G z D (x,y)+ [ t \p a (t,x,y)-e- 2J « t p z (t,x,y))dt + c 5 t- d/2 6 D (x)5 D (y) 



< G z D (x, y) + c 6 t 2 Q - d/2 + est" d/2 5 D (x)5 D (y) 

< G z D (x,y) + c 7 (5 D (x)5 D (y)) 1 - d / 4 . 



(7.11) 
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Since Gp(x, y) > cs{5d{x)8d{v)) 1 d ^ 2 for \x— y\ 2 < SD(x)SD(y) by Theorem ll.il we have from (|7.1ip 
that 

G a D (x,y) < G z (x,y)+c 9 (5 D (x)5 D (y)) d ^G a D (x,y) if \x - y\ 2 < 6 D (x)S D (y). (7.12) 
On the other hand, applying Lemma 17.81 and then Proposition 17.10] we get 

G a D (x,y) < G z D {x,y) + c 10 



S D (x)S D (y) 



\x — y\ d p 

\x-y\< 



Since c\\ ^SAgM^/) < Gjy(x,y) for \x — y\ 2 > <5d(x)<5d(?/) by Theorem ll.lt we have 



G a D (x, y) < G z (x, y) + c 12 \x - y\ p G a D (x, y). (7.13) 
Now using ()7.12p - (l7.13p . we can choose e G (0, 1) small so that 

c 12 G^(x,y)|x-# < ^(x,y) if 8 D (x)8 D (y) < \x - y\ 2 < e 

and 

c 9 (5 D (x)5 D {y)) d ' A G a D (x, y) < \g° d (x, y) if \x - y\ 2 < 5 D (x)5 D (y) < e. 

Thus in these cases, Gp(x,w) < 2Gf)(x,w). 

For the remaining case 5d{x)5d{v) > £, we use (|7.1ip . Lemma 17^31 and (|7.5p to get that 

G a D (x,w) < (l + c 13 (5 D (x)5 D (y)r d / i )G z D (x,y) < (I + c 13 e- d / 4 )G z D (x,y). 

The proof of the theorem is now complete. □ 

We now show that the theorem above covers the case of the independent sum of a Brownian 
motion and a relativistic stable process, and the case of the independent sum of a Brownian motion 
and a truncated stable process. 

For any m > 0, a relativistic a-stable process Y m in M. d with mass m is a Levy process with 
characteristic function given by 



exp( -t ( (|e| 2 +m 2 / Q ) Q/2 -m) ) , £ G ™ 



Suppose that Y m is independent of the Brownian motion X°. We define Z m by := X 4 ° + Y t m . 
We will call the process Z m the independent sum of the Brownian motion X° and the relativistic 
a-stable process Y m with mass m. The Levy measure of Z m has a density 

J zm (x) =A(d, -a)\x\- d - a ^{m x ^\x\) 

where 

m := 2 -(^) r 1 J™ a *F-i e -i-£ ds , 

which is a decreasing smooth function of r 2 . (see |15} pp. 276-277] for details). Thus 

0<jHx)-j zm (x) <c\x\ 2 ~ a ~ d . 

Moreover, the conditions (Al), (A2) and (A3) (a) can be checked easily. Therefore as a corollary of 
Theorem 17. 131 we have the following. 
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Corollary 7.14 There exists a constant C45 = C^(D,a) > such that 

C^GUx,y) < G z D m (x,y) < C 45 G^(x,y), x,y G D, 
where Gf^(x,y) is the Green function of Z m in D. 

By a A-truncated symmetric a-stable process in M. d we mean a pure jump symmetric Levy 
process Y x = (Y t x ,t > 0,P x ,x G R d ) in R d with the Levy density A(d, -a)\x\- d ~ a 1{|x|<a}- Note 
that the Levy exponent tp x of Y x , defined by 



E, 

is given by 



e -ti> (0 f or every x G M d and £ G R d , 



^(0 = -a) / 1 (7-14) 

i{| y |<A} \y\ 

Suppose that Y x is a A-truncated symmetric a-stable process in R d which is independent of the 
Brownian motion X°. We define X x := X® + Y x for i > 0. Then X x has the same distribution as 
the Levy process obtained from X 1 by removing jumps of size larger than A. 

Unlike the symmetric stable process Y, the process Y x can only make jumps of size less than A. 
In order to guarantee the strict positivity of the transition density p^ (t,x,y) for X X ' D , we need 
to impose the following assumption on D. 

Definition 7.15 We say that an open set D in M. d is X-roughly connected if for every x,y G D, 
there exist finitely many distinct connected components Ui, ■ ■ ■ , U m of D such that x G U\, y G U m 
and dist(Uk, E/fe+i) < A for 1 < k < m — 1. 

The following result is proved in [24j . 

Proposition 7.16 (|24 s Proposition 4.4]) For any bounded \-roughly connected open set D in 
M. d , the transition density p^ X (t, x,y) for X X,D is strictly positive in (0, 00) x D x D. 

The other conditions (Al), (A2) and (A3)(b) can be checked easily. Therefore as a corollary of 
Theorem 17. 131 we have the following. 

Corollary 7.17 Suppose D is a bounded X-roughly connected C 1 ' 1 open set in M. d , d > 1. There 
exists C46 = C4Q(D,a) > such that 

C^G l D {x,y) < Gf(x,y) < C 46 G^(x,y), x,y G D, 

where Gp X (x,y) is the Green function of X x in D. 
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